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ABSTRACT 

We present numerical solutions of the 2D relativistic hydrodynamics equations describing the decel- 
eration and expansion of highly relativistic conical jets, of opening angles 0.05 < Oq < 0.2, propagating 
into a medium of uniform density. Jet evolution is followed from a collimated relativistic outflow 
through to the quasi-spherical non-relativistic phase. We show that relativistic sideways expansion 
becomes significant beyond the radius re at which the expansion Lorentz factor drops to 0^^. This 
is consistent with simple analytic estimates, which predict faster sideways expansion than has been 
claimed based on earlier numerical modeling. For t > ts = re / c the emission of radiation from the 
jet blast wave is similar to that of a spherical blast wave carrying the same energy. Thus, the total 
(calorimetric) energy of GRB blast waves may be estimated with only a small fractional error based 
out > ts observations. 

Subject headings: gamma rays: bursts -hydrodynamics -methods: numerical - relativity - radiation 
mechanisms: nonthermal 



1. INTRODUCTION 

The dynamics of gamma-ray burst (GRB) jets and the 
spectral and temporal evolution of their afterglows re- 
main an important problem (see Granot & Ramirez-Ruiz 
2010 for a review). Much of what we know about GRB 
progenitors and the central engines that power them 
comes from multi- wavelength observations of their after- 
glows. From the analytic models which predict the evolu- 
tion of these light curves it is possible to extract estimates 
of jet opening angles, the energetics of the outfiows, and 
the properties of the circumburst medium (Panaitescu & 
Kumar 2002; Cenko et al. 2010; Yost et al. 2003). How- 
ever, more recent numerical modeling has claimed that 
there are strong discrepancies between the analytic and 
numerical models of GRB jets. 

A conical jet-like outfiow expanding at a Lorentz fac- 
tor r evolves as if it were a conical section of a spherical 
outfiow as long as F > 6>^"^, since for V > 6'^^ the (rest 
frame) transverse light crossing time of the jet is larger 
than the expansion/deceleration time. At this stage, the 
fiow is described by the spherical Blandford-McKee (BM) 
blast wave solutions (Blandford & McKee 1976), with 
p2^3 _ \'j Ei^o/ 1^71 nm pC^ . Here r is the blast wave ra- 
dius, n is the ambient medium number density, and £^iso 
is the isotropic equivalent energy, related to the (two- 
sided) jet energy by £^jet = (note that Oq is the 
angular radius). For T > 6^^ a distant on axis observer 
cannot distinguish a jet from a sphere since the emitted 
radiation is beamed into a l/F cone. 

The fiow decelerates to F = at source frame time 



te = re/c = 2?>^ (^iso,53/no)'/' ^, 



^0,-1 day. 
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corresponding to an observer's frame time 

te.^ = -^,te = 0.6 <^ day, (2) 

where Oq = 10-^6>o -i, ^iso = 10^^£^iso,53 erg and n = 
Inocm"^. (For a burst located at redshift all observed 
times should be increased by a factor 1 + z; we do not 
explicitly show this correction in our eqs.) The side- 
ways expansion is expected to be relativistic as long as 
the blast wave is relativistic and the post-shock energy 
density is relativistic (Rhoads 1999). If this is the case, 
ai t > te the lateral expansion rapidly increases the jet 
opening angle and accelerates its deceleration (Rhoads 
1999), reducing F to ~ 1 with only a logarithmic in- 
crease of r (to ~ In^^"*^ X re). Thus, the observed time 
scale for the fiow to become transrelativistic is (Livio & 
Waxman 2000) 



ts,(B ^ te,(B -^re/c^te = 230 (£^iso, 53/^0) 0\ 



1/3 ^2/3 
0, 



-1 day. 
(3) 

On a similar time scale, the fiow is expected to become 
quasi-spherical, i.e. the jet is expected to expand to ^ ~ 
l,and the outfiow is subsequently expected to evolve into 
the spherical non-relativistic Sedov-von Neumann- Taylor 
(ST) fiow. 

This simple analytic description of jet expansion was 
challenged by a series of numerical calculations (Granot 
et al. 2001; Cannizzo et al. 2004; Zhang & MacFadyen 
2009; Meliani & Keppens 2010). It was argued, based on 
the numerical results, that the sideways expansion of the 
jet is not relativistic, and that the jet retains its narrow 
original opening angle, 6>o, as long as it is relativistic 
(Granot 2007; Zhang & MacFadyen 2009). This implies 
that the jet continues to evolve like a conical section of 
a spherical outfiow with energy E\^o^ with Lorentz factor 
following the BM solution, up to the radius tnr = c^nr 
at which it becomes sub-relativistic, 
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The different descriptions of jet evolution inferred from 
analytic and numeric modeling lead to different predic- 
tions for the observed properties of GRB afterglows (e.g. 
van Eerten et al. 2010a). For example, the suppression 
of the observed flux produced by the jet blast wave at 
t > t6',e, compared to that produced by a spherical blast 
wave with the same E^so^ is smaller (i.e. the "jet break" 
is less pronounced) if the jet does not expand signifi- 
cantly while it is relativistic. Furthermore, if the jet 
does not expand significantly and remains highly colli- 
mated and relativistic at t > , then the accuracy of the 
late-time calorimetric estimates of the jet energy, which 
assume quasi-spherical emission at t ~ (Frail et al. 
2000; Berger et al. 2004), is questionable. 

The main goal of this letter is to resolve the apparent 
discrepancy between the analytic and numeric descrip- 
tion of jet expansion. A more detailed discussion of the 
properties of the flow and of the emitted radiation will 
be given in a subsequent more comprehensive publication 
(Wygoda & Waxman 2011). Our numerical calculations 
are described in § 2, and their results regarding jet ex- 
pansion are described in § 3. In § 4 we briefly discuss the 
implications for jet breaks and GRB calorimetry. Our 
conclusions are summarized in § 5. 

2. NUMERICAL CALCULATIONS 

We use the RELDAFNA code (Klein 2010) to nu- 
merically solve the 2D special relativistic hydrodynam- 
ics equations describing the flow of an ideal fluid with 
a constant poly tropic index, 7 = 4/3. RELDAFNA is 
a Godunov type Eulerian code, with second order accu- 
racy time and space integration. It uses adaptive mesh 
refinement and is massively parallelized, allowing the use 
of effectively high resolution even in multiscale problems 
such as the current jet simulations. RELDAFNA was 
tested (Klein 2010, see also § 3) by comparing its solu- 
tions for standard test problems to those of similar codes 
(Zhang & MacFadyen 2006; Meliani et al. 2007), and was 
shown to perform similarly. 

The initial conditions chosen for our numerical calcu- 
lations were a conical section of opening angle ^0 within 
which the flow fields are given by the BM solution for 
^iso = 10^^ erg and initial density n = lcm~^, sur- 
rounded by a static uniform cold gas of density n = 
lcm~^ and pressure po = lO~^^nmpC^ . The radius 
of the conical section was chosen so that the Lorentz 
factor of the fluid behind the shock is F = 20. We 
present solutions for ^0 =0.2, 0.1 and 0.05 (correspond- 
ing to E^et = 2 • 10^\5 • 10^^ and 1.25 • 10^° erg). The 
6^0 = 0.2 simulation is similar to the simulation presented 
in (Zhang & MacFadyen 2009, hereafter ZM09). The 
only difference is that we use 7 = 4/3, instead of an 
equation of state for which 7 varies smoothly between 
7 = 4/3 for relativistic material and 7 = 5/3 for non- 
relativistic material. Our choice is inaccurate for late 
times, when the flow becomes non-relativistic, but this 
inaccuracy is not expected to affect our results qualita- 
tively. Moreover, if a significant fraction of the post- 
shock energy density is carried by magnetic fields and 
relativistic electrons, as required in order to account for 
afterglow observations, 7 remains close to its relativistic 
value further into the non-relativistic flow stage. 

The size of the finest numerical cells in the simulation 
was initially taken as 5.6 • 10^^ cm, similar to ZM09. The 



results of our simulations were checked for convergence 
by increasing the grid resolution by a factor of 4 in each 
dimension. Increased resolution calculations were car- 
ried out both for initial conditions identical to those of 
the nominal calculations, and for initial conditions with 
a reduced radius of the conical section corresponding to 
a post-shock BM Lorentz factor of F = 40. These con- 
vergence tests indicated that while the Lorentz factor 
behind the shock, as well as the lightcurves of the high 
frequencies that depend strongly on the high F region, 
are not yet fully converged, the spreading of the jet is 
converged to a level of 10%. For example, the time it 
takes the jet angle to double its initial value decreased 
in the convergence test by ~ 8% for ^0 = 0.2 and by 
~ 25% for Oq = 0.05. The results presented in the next 
sections are the ones obtained with the higher numerical 
resolution. A more detailed analysis of the numerics will 
be given in Wygoda & Waxman (2011). 

3. JET EXPANSION 
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Fig. 1. — The evolution of the jet opening angle, ^(t)/^o, as a 
function of time (measured in the source frame and normalized to 
te), for ^0 = 0.2 (solid), 0.1 (dashed) and 0.05 (dash-dotted). is 
defined as the cone opening angle within which 95% of the energy, 
excluding rest mass energy, is included. The 3 lines denoted tNR 
show ^nr/^^ for the 3 values of ^o- 

Figure 1 presents the evolution of the jet opening angle, 
0(t\ as a function of time, d is defined as the cone open- 
ing angle within which 95% of the energy, excluding rest 
mass energy, is included. We find that significant side- 
ways expansion begins at t ~ ^51 <C ^nr, in accordance 
with the analytic estimates described in § 1. Narrower 
jets begin expanding earlier, Iq oc in accordance 

with the analytic estimates, and are therefore expected 
to also decelerate earlier. 

The latter point is demonstrated in figure 2, which 
shows the density distribution of the ^0 =0.2 and 0.05 jet 
flows at identical time t = 0.95tNR (note that ^nr is inde- 
pendent of 6>o). The 6^0 =0.2 jet has tripled its opening 
angle and its tip is still close to its "isotropic equivalent 
location" , i.e. the location of a spherical blast wave with 
the same ^iso- The opening angle of the 6>o = 0.05 jet 
has increased by more than an order of magnitude, and 
its spreading has significantly slowed down its expansion. 
The influence of jet expansion at ^nr is much stronger for 
the narrower jet, in accordance with the analytic anal- 
ysis described in § 1: The ratio of ^nr to tg is close to 



Jet Dynamics and Calorimetry 



3 




-1 -0.5 0.5 1 
x[10^^cm] 



Fig. 2. — The density distribution of the Oq = 0.2 (right) and 
0.05 (left) jet flows at identical time t = 0.95tNR (^nr is the same 
for both jets since Eiso and n are identical for both). The black 
lines indicate the shape and size the jet would have reached if it 
behaved like a conical section of a sphere. 

unity for the = 0.2 jet, and significantly larger for the 
= 0.05 jet: at the source frame tg ^ to -\- tq/c = 2to so 

that ^nrA, ^ l/20l^^ = 1.5 ,3.7 for Oq = 0.2 ,0.05. 
We thus find that the jet sideways expansion is rela- 

2 /3 

tivistic and becomes significant dit t ^ to (x Oq , and that 
this expansion leads to deceleration to sub-relativistic ve- 

2 /3 

locity on a time scale tg (xOq , which for <^ 1 is much 
smaller than ^nr c<: 6q. This behavior is consistent with 
the analytic analysis, and inconsistent with the claims 
based on earlier numerical modeling, that jet expansion 
is not significant up to t ~ ^nr- 
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Fig. 3. — Comparison of Q as function of source frame time ob- 
tained for the ^0 = 0.2 jet in our numerical calculation (solid line), 
with that obtained by ZM09 (dotted), and by the modified Rhoads' 
model (dashed, eq. 5). The short vertical line denotes the time at 
which r = 2, up to which Rhoads' model applies. The dash-dotted 
line shows the exponential, In^/^o ^ t: to which the numerical 
results were compared in ZM09 (see their fig. 3). 

In order to identify the origin of this apparent discrep- 
ancy, we compare our numerical results to those of ZM09 



in fig. 3. The figure demonstrates that the jet expansion 
obtained in our calculation is similar to that obtained in 
ZM09. The conclusion that sideways expansion is not 
relativistic, and unimportant until t ^ ^nr, was reached 
in ZM09 based on noting that the growth of 0{t) is much 
slower than the exponential growth expected for relativis- 
tic sideways expansion (see fig. 3). This conclusion is, 
however, not valid, since exponential growth is expected 
only for ^ r ^ 1, and is not applicable for the evo- 
lution of the 6>o = 0.2 jet under consideration, for which 
expansion becomes significant only for F < = 5. 

For this regime of F, one cannot use the exponential 
approximation, but should rather solve the differential 
equation 

de/dr = Cs/Tcr, (5) 

describing relativistic sideways expansion at the post 
shock speed of sound in the jet frame (c^ = c/>/3 for 
F 1), along with mass and momentum conservation, 
that determine F(r, 0) (for more details see Rhoads 1999). 
Note that we replace eq. (3) of Rhoads (1999) with the 
more accurate eq. (5) (see also Piran 2000). This mod- 
ification leads to a significant modification of 0{r) only 
for > 0.4. For (9o"^ > F > 1, the solution of eq. (5) 
is indeed exponential, IuO/Oq oc r^/^. However, such a 
regime does not exist for 6>o = 0.2. As demonstrated in 
the figure, the numerical results are in good agreement 
with the solution of the simple model of eq. (5), and 
therefore confirm the validity of the analytic estimates 
described in § 1. 

4. LIGHT CURVES AND CALORIMETRY 

We have calculated the synchrotron emission expected 
to be produced by shock accelerated electrons assuming 
that the magnetic field and the electrons hold a constant 
fraction eg = = 0.1 of the internal energy, and that 
the electron energy distribution is a power law with index 
p = 2.4. Electron cooling and synchrotron self absorption 
are neglected. 

Figure 4 shows radio light-curves {u ~ 3GHz, for which 
self-absorption is not important) predicted by the nu- 
merical model for the = 0.2 jet, for observers lying on 
the jet axis and at an angle ^obs. = ^o- The numerical 
light curves of the jet are compared with those predicted 
for spherical (ID) fireballs with total energy £^iso and 
Ejet, as well as with that predicted for a conical section 
with opening angle of a spherical E[so fireball, repre- 
senting the lightcurve predicted for a non-expanding jet. 
Also shown is the radio fiux of a spherical fireball with 
energy E'jet, assuming its evolution is described by the 
non-relativistic Sedov-von Neumann- Taylor solution. 

The jet emission is suppressed, compared to that of a 
spherical blast wave with energy Ei^o^ dit t > to. The 
suppression is larger than would be predicted for a non- 
expanding jet (i.e. not due to the "missing fiux" from 
the absent > Oq parts of the shell, but rather to the jet 
spreading, in accordance with van Eerten et al. 2010b). 
The figure also demonstrates that at t > , the observed 
fiux is similar to that of a spherical fireball with energy 
Ejet, for observers lying both on- axis and at an angle 
= Oq. Moreover, although at t ~ the jet has not yet 
reached full spherical symmetry and is still mildly colli- 
mated (see fig. 1), at t > ts the fiux is well approximated 



4 



Wygoda et al. 



by that of a blastwave following the non-relativistic ST 
evolution. 




10 10 

Fig. 4. The (2D) Oq = 0.2 jet radio lightcurve {y ^ GHz) 
for observers lying on-axis (solid black line) and at 0^ = 0.2 rad, 
compared with those of several spherical (ID) fireballs: E = E[s 



(dashed), E ■ 
ted), and E - 
totic ST behavior. 



E[sQ with emission only from conical section (dot- 
E'jet (dash-dotted). The red line denotes the asymp- 



In Figure 5 we compare the Oq = 0.2,0.1,0.05 jet 
lightcurves to those of spherical (ID) fireballs with the 
corresponding E = £^jet- Scaling with ts,^ brings the 
lightcurves to a similar (universal) form, implying that 
significant jet spreading occurs dit t ^ te (x d^l'^ . Al- 
though, as noted above, at t ~ the jet does not yet 
reach full spherical symmetry, from ~ O.Stg onward the 
lightcurves do not depart from those of the ID = Ejet 
fireballs by more than ~ 50%. A more detailed analy- 
sis of the jet lightcurves will be presented in Wygoda & 
Waxman (2011). 

5. DISCUSSION AND CONCLUSIONS 

While there still remain many unsolved problems re- 
garding the structure and dynamics of GRB jets that can 
only be addressed by detailed hydrodynamic modeling, 
our work has demonstrated that analytic estimates pro- 
vide a reasonable description of the behavior of the jet 
and the evolution of its afterglow. 

We have shown that relativistic sideways expansion be- 
comes significant at t > tgi ex 6>q^^ (fig. 1), in accordance 
with analytic estimates, and that the expansion is well 
described by the modified Rhoads model (fig. 3, eq. 5)). 
Our numerical results are consistent with those of ZM09, 
who calculated = 0.2 jet evolution. The apparent dis- 
crepancy between earlier numerical and analytic results 



arose because the simulations weren't compared to the 
full solution of Rhoads' model (fig. 3, § 3), and because 
for the large 6>o chosen it is difficult to test the relativis- 
tic expansion assumption since ^nr and is are similar, 

^nrA. ^ 1/2^^/' = 1.5 (see § 3). 

Jet expansion has a significant effect on its observed 
properties. The suppression of the ffux at t > ^6* is 
stronger than in the absence of spreading (fig. 4), and 
at t > ^5 the emission of radiation from the jet blast 
wave is similar to that of a spherical blast wave carry- 
ing the same energy (fig. 4,5). Moreover, although at 
t ^ is the jet has not yet reached full spherical symme- 
try and is still mildly collimated (see fig. 1), at t > tg 
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Fig. 5. — Ratio between the 2D jet flux and that of a ID spherical 
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the fiux is well approximated by that of a blastwave fol- 
lowing the non-relativistic ST evolution (fig. 4). Thus, 
the total (calorimetric) energy of GRB blast waves may 
be estimated with only a small fractional error based on 
t>ts observations. We expect to see this technique and 
its variants (Frail et al. 2000; van der Horst et al. 2008; 
Shivvers & Berger 2011) applied to increasing numbers 
of GRB afterglows when the new generation of facilities 
(EVLA, LOFAR) starts fuh operation. 
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